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Hypergraph states form a family of multiparticle quantum states that generalizes the well-known 
concept of Greenberger-Horne-Zeilinger states, cluster states, and more broadly graph states. We 
study the nonlocal properties of quantum hypergraph states. We demonstrate that the correlations 
in hypergraph states can be used to derive various types of nonlocality proofs, including Hardy- 
type arguments and Bell inequalities for genuine multiparticle nonlocality. Moreover, we show that 
hypergraph states allow for an exponentially increasing violation of local realism which is robust 
against loss of particles. Our results suggest that certain classes of hypergraph states are novel 
resources for quantum metrology and measurement-based quantum computation. 
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Introduction .— Multiparticle entanglement is central 
for discussions about the foundations of quantum me¬ 
chanics, protocols in quantum information processing, 
and experiments in quantum optics. Its characterization 
has, however, turned out to be difficult. One problem 
hindering the exploration of multiparticle entanglement 
is the exponentially increasing dimension of the Hilbert 
space. This implies that making statements about gen¬ 
eral quantum states is difficult. So, one has to concen¬ 
trate on families of multiparticle states with an easier-to- 
handle description. In fact, symmetries and other kinds 
of simplifications seem to be essential for a state to be 
a useful resource. Random states can often be shown to 
be highly entangled, but useless for quantum information 
processing [1]. 

An outstanding class of useful multiparticle quantum 
states is given by the family of graph states [2], which 
includes the Greenberger-Horne-Zeilinger (GHZ) states 
and the cluster states as prominent examples. Physically, 
these states have turned out to be relevant resources 
for quantum metrology, quantum error correction, or 
measurement-based quantum computation [2]. Math¬ 
ematically, these states are elegantly given by graphs, 
which describe the correlations and also a possible inter¬ 
action structure leading to the graph state. In addition, 
graph states can be defined via a so-called stabilizer for¬ 
malism: A graph state is the unique eigenstate of a set of 
commuting observables, which are local in the sense that 
they are tensor products of Pauli measurements. These 
stabilizer observables are important for easily comput¬ 
ing correlations leading to violations of Bell inequalities 
[3, 4], as well as designing simple schemes to characterize 
graph states experimentally [5]. 

Recently, this family of states has been generalized to 
hypergraph states [6-11]. These states have been rec¬ 
ognized as special cases of the so-called locally maxi¬ 
mally entangleable (LME) states [6]. Mathematically, 
they are described by hypergraphs, a generalization of 
graphs, where a single hyperedge can connect more than 
two vertices. They can also be described by a stabilizer 





FIG. 1: Examples of hypergraphs, (a) A simple hypergraph 
with three vertices and a single edge e = {1, 2, 3} connecting 
all three vertices. The corresponding hypergraph state IHs) = 
Gi 23 |-I-)®® = (|000)-l-|001}-l-- • --l-lllO} —1111))/-^^^ discussed 
in detail in the text, (b) This hypergraph contains only two- 
edges, so it is an ordinary graph. The state corresponding to 
this fully connected graph is the five-qubit GHZ state, (c) 
The fully connected three-uniform hypergraph represents a 
state that can be seen as a generalization of the GHZ state. 


formalism, but this time, the stabilizing operators are 
not local. So far, hypergraph states have turned out to 
play a role for search algorithms in quantum comput¬ 
ing [12], quantum fingerprinting protocols [13], and they 
have been shown to be complex enough to serve as wit¬ 
nesses in all QMA problems [14]. They have recently 
been investigated in condensed matter physics as ground 
states of spin models with interesting topological proper¬ 
ties [15, 16]. In addition, equivalence classes and further 
entanglement properties of hypergraph states have been 
studied [9]. 

In this paper we show that hypergraph states violate 
local realism in an extreme manner, but in a way that is 
robust against loss of particles. We demonstrate that this 
leads to applications of these states in quantum metrol¬ 
ogy and quantum computation. We see that the sta¬ 
bilizer formalism describing hypergraph states, despite 
being nonlocal, can be used to derive Hardy-type nonlo¬ 
cality arguments [17], Bell inequalities for genuine multi¬ 
particle entanglement [18], or a violation of local realism 
with a strength exponentially increasing with the num¬ 
ber of particles. Our approach starts precisely with the 
properties of the stabilizer, in order to identify the usefnl 
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correlations provided by quantum mechanics. This is in 
contrast to previous approaches that were either too gen¬ 
eral, e.g. Bell inequalities for general multiparticle states 
[19, 20], or too restricted, considering only few specific 
examples of hypergraph states and leading to non robust 
criteria [9]. The violation of local realism is the key to 
further applications in information processing: Indeed, it 
is well known that violation of a Bell inequality leads to 
advantages, in distributed computation scenarios [ 21 , 22 ]. 
In addition, we will explicitly show that certain classes of 
hypergraph states lead to Heisenberg scaling in quantum 
metrology and advantages in measurement-based quan¬ 
tum computation. 

Hypergraph states .— A hypergraph H = {V,E) con¬ 
sists of a set of vertices V = {1,..., N} and a set of hy¬ 
peredges E C 2^, with 2^ the power set of V. While for 
graphs edges connect only two vertices, hyperedges can 
connect more than two vertices; examples of hypergraphs 
are depicted in Fig. 1. For any hypergraph we define the 
corresponding hypergraph state \H) as the Wqubit state 

\H) = n C'e|+)®^, (1) 

eeE 

where J-l-) = (JO) -I- \l))/\/2, e is a hyperedge and Ce 
is a multi-qubit phase gate acting on the Hilbert space 
associated with the vertices v G e, given by the matrix 
Ce = 1 — 2jl... 1)(1... Ij. The first nontrivial hyper¬ 
graph state consists of = 3 qubits connected by a 
single hyperedge [see Fig. 1(a)]. Hypergraph states have 
been recognized as special cases of LME states, generated 
via a fixed interaction phase of (/> = tt [ 6 ] . 

Alternatively, we can define the hypergraph states us¬ 
ing a stabilizer formalism [9]. For each qubit i we define 
the operator 

gi = Xi^Ce\{i}. ( 2 ) 

eeE 

Here and in what follows, we denote by Xi and Zi the 
Pauli matrices, acting on qubit. The hypergraph state 
can be defined as the unique eigenstate for all of them, 
gi\H) = \H) with the eigenvalue -1-1. Consequently, the 
hypergraph state is an eigenstate of the entire stabilizer, 
i.e., the commutative group formed by all the products 
of the gi- It should be noted that the gi are, in general, 
non-local operators, as they are not tensor-products of 
operators acting on single parties. We say that a hyper¬ 
edge has cardinality fc, if it circumscribes k vertices and a 
hypergraph is k-uniform, if all edges are fc-edges. Finally, 
note that different hypergraphs may lead to equivalent 
hypergraph states, in the sense that the two states can 
be converted into one other by a local basis change. For 
small numbers of qubits, the resulting equivalence classes 
have been identified [9]. 

Loeal correlations from the nonlocal stabilizer .— The 
key observation for the construction of our nonlocality 


arguments is that the stabilizer of hypergraph states, 
despite being nonlocal, predicts perfect correlations for 
some local measurements. In the following, we explain 
this for the three-qubit hypergraph state li^ 3 ), but the 
method is general. The stabilizing operators for the 
three-qubit hypergraph state are 

3i = All 0 C 23 , 52 = X 2 0 (713, ^3 = 0 C 12 . (3) 

We can expand the controlled phase gate Cy on two 
qubits, leading to 

51 = 0 (100)(001 + 101)(011 + 110)(101 - 111)(111) (4) 

and similar expressions for the other gi. Since gi\Hf) = 
+ \H 3 ), the outcomes for X measurements on the first 
qubit and Z measurements on the second and third 
qubits are correlated: if one measures ” -|- ” on the first 
qubit, then the other two parties cannot both measure 
” — ” in Z direction, as this would produce —1 as the 
overall eigenvalue. So, we extract the first correlation 
from the stabilizer formalism: 

P{+--\XZZ) = 0. (5) 

The l.h.s of Eq. (5) denotes the probability of measuring 
-I-in XZZ on the qubit 1, 2, and 3, respectively. Sim¬ 

ilarly, it follows that if one measures ” — ” in Af direction 
on the first qubit, then the other parties, both have to 
measure ” — ” in Z direction. So we have: 

P{- + + \XZZ) + P{- + -\XZZ) + P{- -+\XZZ) = 0, 

_ ( 6 ) 

which implies, of course, that each of the probabilities 
is zero. Since the three-qubit hypergraph state is sym¬ 
metric, the same correlations for measuring X on other 
qubits can be obtained by considering (72 and gs, leading 
to permutations of correlations in Eq. (5 , 6 ). 

The three-qubit hypergraph state {H ^).— We start with 
the discussion of fully local hidden variable (HV) models. 
Such models assign for any value of the HV A results to all 
measurements of the parties in a local manner, meaning 
that the probabilities for a given HV factorize. If we 
denote by the result and by Si the measurement setting 
on the particle, respectively, then the probabilities 
coming from local models are of the form 

H(ri,r2,r3lsi,S2,S3) = (7) 

= J dXp{X)x^{ri\si, X)x^{r 2 \s 2 , X)X^{r3\s3, X). 

Eor probabilities of this form, it is well known that it 
suffices to consider models which are, for a given A, de¬ 
terministic. This means that takes only the values 0 
or I, and there is only a finite set of y® to consider. 
Observation 1. If a fully local hidden variable model 
satisfies the conditions from Eq. (5, 6) and their sym¬ 
metric correlations coming from the permutations, then 
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it must fulfill 

P(+- \XXX)+P{-+-\XXX)+P{ - + \XXX) = 0. 

( 8 ) 

The proof of this statement is done by exhausting all 
possible local deterministic assignments. 

In contrast, for \Hf) we have 

P[+--\XXX) = ^ (9) 

and the same holds for the permutations of the qubits. 
The above is a so-called Hardy argument [17], namely, 
a set of joint probabilities equal to 0 or, equivalently, 
logical implications that together implies that some other 
probability is equal to zero. 

Our method shows how the correlations of the nonlo¬ 
cal stabilizer can be used for Hardy-type arguments. We 
recall that Hardy-type arguments have been obtained for 
all permutation-symmetric states [23, 24]. However, they 
involved different settings and have no direct connection 
with the stabilizer formalism, making a generalization 
complicated. In contrast, we will see that our measure¬ 
ments can even be used to prove genuine multiparticle 
nonlocality of the hypergraph state. First, we translate 
the Hardy-type argument into a Bell inequality: 
Remark 2. Putting together all the null terms derived 
from the stabilizer formalism and subtracting the terms 
causing a Hardy-type argument, we obtain the Bell in¬ 
equality 

= [P(+ - -\XZZ) + P{- + +\XZZ) 

-\-P {—I— \XZZ) -\- P{ - hlXZZ)-!- permutat.J 

— [P(-|- \XXX)-\- permutations] >0, (10) 

where the permutations include all distinct terms that are 
obtained by permuting the qubits. The three-uniform hy¬ 
pergraph state violates the inequality (10) with the value 

of (4'^) = -3/16. 

This Bell inequality follows from the Hardy argument: 
If a deterministic local model predicts one of the results 
with the minus signs, it also has to predict at least one of 
the results corresponding to the terms with a plus sign, 
otherwise it contradicts with the Hardy argument. In 
addition, all the terms with a minus sign are exclusive, 
so a deterministic LHV model can predict only one of 
them. 

The Hardy-type argument and the Bell inequality can 
be generalized to a higher number of qubits, if we con¬ 
sider fV-qubit hypergraphs with the single hyperedge 
having a cardinality N: 

Observation 3. Consider the N-qubit hypergraph state 
with a single hyperedge of cardinality N. Then, all the 
correlations coming from the stabilizer [as generalizations 
of Eqs. (5,6)] imply that for any possible set of results 
{ri} where one ri.^ = -|-1 and two ri^ = = — 1 one has 

P(ri,r2,...,r^lXiX2...Xv) = 0. (11) 


For the hypergraph state, however, this probability equals 
This Hardy-type argument leads to a Bell in¬ 
equality as in Eq. (10) which is violated by the state with 
a value of-{2^ - N - 2). 

Clearly, the violation of the Bell inequality is not 
strong, as it does not increase with the number of parti¬ 
cles. Nevertheless, Observation 3 shows that the nonlocal 
stabilizer formalism allows one to easily obtain nonlocal¬ 
ity proofs. In fact, one can directly derive similar argu¬ 
ments for other hypergraph states (e.g. states with one 
hyperedge of cardinality N and one further arbitrary hy¬ 
peredge), these results will be presented elsewhere. Note 
that these states are not symmetric, so the results of 
Refs. [23, 24] do not apply. 

So far, we only considered fully local models, where 
for a given HV all the probabilities factorise. Now 
we go beyond this restricted type of models to the 
so-called hybrid models [18]. We consider a bipar¬ 
tition of the three particles, say A\BC, and con¬ 
sider a model of the type P(ri, r 2 ,rajsi, S 2 , S 3 ) = 

/dAp(A)x^(rilsi,A)x®‘^(r 2 ,r 3 ls 2 ,S 3 , A). Here, Alice is 
separated from the rest, but may contain correla¬ 
tions, e.g., coming from an entangled state between B 
and C. In order to be physically reasonable, however, we 
still request x^*^ not to allow instantaneous signaling. 

This kind of models, even if different bipartitions are 
mixed, cannot explain the correlations of the hyper¬ 
graph state, meaning that the hypergraph state is gen¬ 
uine multiparticle nonlocal. First, one can see by direct 
inspection that the stabilizer conditions from Eqs. (5, 
6 ) are not compatible with the hypergraph correlations 
P(- \XXX) = 1/16 and P(- \ZZZ) = 1/8. Con¬ 

trary to the correlations in Eq. (9) these are symmetric, 
and allow the construction of a Bell-Svetlichny inequality 
[18] valid for all the different bipartitions: 

Observation 4. Putting all the terms from the hyper¬ 
graph stabilizer formalism and the correlations P(- 

— \XXX) and P(- \ZZZ) together, we obtain the 

following Bell-Svetlichny inequality for genuine multipar¬ 
ticle nonlocality, 

= [P(+ - -\XZZ) + P(- + +\XZZ) 

-|-P(—I— \XZZ)-\-P{ - \-\XZZ)-\- permutat.] 

+P{ - \XXX)-P{ - \ZZZ)>0, (12) 

which is violated by the state jPs) with {B^ 7 = —1\16. 

The proof is done by an exhaustive assignments of 
nonsignaling and local models. 

To investigate the noise tolerance of Ineq. (12), we con¬ 
sider states of the type g = (1 — e)lP)(Pl-l-el/8 and ask 
how much noise can be added, while the inequality is still 
violated. The white noise tolerance of the Ineq. (12) is 
£ = 1/13 ~ 7.69% and is optimal in the sense that for 
larger values of e a hybrid model can be found which 
explains all possible measurements of X and Z (within 
numerical precision). The existence of such a model can 
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be shown by linear programming (see Appendix A [25]). 
With the same method we can also prove that the state 
becomes fully local with respect to X and Z measure¬ 
ments for e > 2/3 ~ 66.6%. 

Three uniform hypergraph states. — Let us extend our 
analysis to hypergraph states with a larger number of 
particles. Here, it is interesting to ask whether the viola¬ 
tion of Bell inequalities increases exponentially with the 
number of parties. Such a behaviour has previously been 
observed only for GHZ states [3] and some cluster states 

[4]. 

GHZ states are described by fully connected graphs 
(see Fig. 1), i.e., fully connected two-uniform hypergraph 
states. It is thus natural to start with fully connected 
three-uniform hypergraph states. First, we observe that 
for such states on N qubits and for even m with 1 < m < 
N 


{X_^Z ...Z) 

m 


+ 4 if TO = 2 mod 4, 
— 4 if TO = 0 mod 4. 


(13) 


Moreover, if to = A^, then the correlations are given by 


( AA...AA ) 

N 


0 if A = 0 mod 4, 
1 if A = 2 mod 4. 


(14) 


Finally, we always have {ZZ... ZZ) = 0 (see Appendix 
B for details [25]). 

We then consider the following Bell operator 

Bn = — [AAA... AA] -|- [HHA ... A + permutat.] — 

— [BBBBA ... a -I- permutat.] -I- [. ■. ] — ... (15) 

Note that this Bell operator is similar to (Hjjf) of the 
original Mermin inequality [3], but it differs in the num¬ 
ber of B (always even) considered. Using the correlations 
computed above, we can state: 

Observation 5. If we fix in the Bell operator Bn in 
Eg. (15) the measurements to be A = Z and B = X, 
then the N-qubit fully-connected three-uniform hyper¬ 
graph state violates the classical bound, by an amount 
that grows exponentially with number of qubits, namely 

{Bn)c < for local HV models, and 

{Bn)o > — - for the hypergraph state. (16) 

V 2 

The proof is given in Appendix G [25]. 

Four-uniform hypergraph states. — Finally, let us con¬ 
sider four-uniform complete hypergraph states. For 
them, the correlations of measurements as in Eq. (13) 
are not so simple: They are not constant, and depend 
on TO as well as on A. Nevertheless, they can be ex¬ 
plicitly computed, and detailed formulas are given in the 
Appendix D [25]. From these correlations, we can state: 
Observation 6. The N-qubit fully-connected four- 
uniform hypergraph state violates local realism by an 


amount that grows exponentially with number of qubits. 
More precisely, one can find a Mermin-like Bell operator 
Bn such that 

( \N-\ 

{Bn)q N—¥00 + _ 1.20711^ 

{^n)c ^ ^ 2 ^^+" 2^2 + 2 ■ 

A detailed discussion is provided in Appendix E [25]. 

Robustness. — So far, we have shown that three- and 
four-uniform hypergraph states violate local realism com¬ 
parable to GHZ states. A striking difference is, however, 
that the entanglement and Bell inequality violation of 
hypergraph states is robust under particle loss. This is 
in stark contrast to GHZ states, which become fully sep¬ 
arable if a particle is lost. We can state: 

Observation 7. The N-qubit fully-connected four- 
uniform hypergraph state preserves the violation of the 
local realism even after loss of one particle. More pre¬ 
cisely, for A = 8fc -I- 4 we have {Bn-i)q/{Bn)q 
l/(-\/2 + 1). This means that the reduced state shows the 
same exponential scaling of the Bell inequality violation 
as the original state. 

Eor the detailed discussions see Appendix F [25]. 

For three-uniform complete hypergraph states we can 
prove that the reduced states are highly entangled, as 
they violate inequalities testing for separability [26] ex¬ 
ponentially. This violation decreases with the number of 
traced out qubits, but persists even if several qubits are 
lost. This suggests that this class of hypergraph states is 
also more robust than GHZ states, details can be found 
in Appendix G [25]. Despite the structural differences, 
this property resembles of the W state, which is itself less 
entangled but more robust than the GHZ state [31]. In 
addition, this may allow the lower detection efficiency in 
the experiments. 

Discussion and conclusion. — A first application of our 
results is quantum metrology. In the standard scheme 
of quantum metrology one measures an observable Mg 
and tries to determine the parameter 9 which describes a 
phase in some basis [27, 28]. If one takes product states, 
one obtains a signal (Mg) ^ cos {9) on a single parti¬ 
cle, repeating it on A particles allows to determine 9 
with an accuracy 59 ^ 1/-\/A, the so-called standard 
quantum limit. Using an A-qubit GHZ state, however, 
one observes {{Mg)®^) ^ cos {N9) and this phase super¬ 
resolution allows to reach the Heisenberg limit 59 ^ 1/A. 
For a general state g, it has been shown that the visibility 
of the phase super-resolution is given by the expectation 
value of the Mermin-type inequality, V = Tr{BN g)/‘2^~^ 
[28]. So, since the three-uniform hypergraph states vio¬ 
late this inequality with a value {Bn)q ~ the 

visibility is U ~ 1/2, independently of the number of 
particles. This means that these states can be used for 
Heisenberg-limited metrology, and from our results they 
can be expected to have the advantage of being more 
robust to noise and particle losses. 
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The second application of exponential violation of Bell 
inequalities is nonadaptive measurement based quantum 
computation with linear side-processing [NMQC^) [29]. 
NMQC^ is a non-universal model of quantum compu¬ 
tation where linear classical side-processing is combined 
with quantum measurements in a nonadaptive way, i.e., 
the choice of settings is independent of previous out¬ 
comes. In Ref. [29] the authors connect the expecta¬ 
tion value of a full-correlation Bell expression [20] with 
the success probability of computing a Boolean function, 
specified as a function of the inequality coefficients, via 
NMQCq. In particular, the exponential violation of gen¬ 
eralized Svetlichny inequalities [30] (equal to Mermin in¬ 
equalities for even TV ), corresponds to a constant success 
probability Psucc of computing the pairwise AND on N 
bits extracted from a uniform distribution, whereas in the 
classical case Psucc ~ 1/2 decrease exponentially with N. 
As a consequence, the exponential violation of the full- 
correlation Bell expression Bn can be directly related to 
an exponential advantage for computation tasks in the 


NMQC(^ framework. Moreover, in several cases, e.g., 
4-uniform hypergraph states of = 6 mod 8 qubits, 
also the Svetlichny inequality is violated exponentially, 
providing an advantage for computation of the pairwise 
AND discussed in Ref. [29]. 

In summary, we have shown that hypergraph states 
violate local realism in many ways. This suggests that 
they are interesting resources for quantum information 
processing, moreover, this makes the observation of hy¬ 
pergraph states a promising task for experimentalists. In 
our work, we focused only on some classes of hypergraph 
states, but for future research, it would be desirable to 
identify classes of hypergraph states which allow for an 
all-versus-nothing violation of local realism or which are 
strongly genuine multiparticle nonlocal. 

We thank D. Nagaj, F. Steinhoff, M. Wiesniak, and 
B. Yoshida for discussions. This work has been sup¬ 
ported by the EU (Marie Curie CIG 293993/ENFOQI), 
the FQXi Fund (Silicon Valley Community Foundation), 
the DFG and the ERG. 


APPENDIX A. GENUINE MULTIPARTICLE NONLOCALITY WITH LINEAR PROGRAMMING 


Fixed the number of measurement settings and outcomes, probabilities arising from a hybrid local-nonsignalling 
model, as the one described in the main text for the splitting A\BC, form a polytope whose extremal points are given 
by combination of deterministic local assignments for the party A and extremal nonsignalling assignments, i.e., local 
deterministic and PR-boxes [21], for the parties BC. In order to detect genuine multiparticle nonlocality, one has 
to consider all combinations of probabilities arising from the other local-nonsignalling splitting, namely C\AB and 
B\AC. Geometrically, this corresponds to take the convex hull of the three polytopes associated with the three different 
splitting. Let us denote such polytopes as 'Pa\bc^'Pc\ab^'Pb\ac S'lid their convex hull as Vl-ns- By definition of 
convex hull, every vector p G Vl-ns can be written as a convex combination of three vectors Pa\bCiPc\ab and 
Pb\aCi which, in turn, can be written as a convex combination of the vertices of the corresponding polytope. 

To check whether a given point p belongs to Vl-ns it is, therefore, sufficient the description in terms of the extremal 
points of Va\bCjVc\ab,Vb\ac- Let us denote them as {v^}. The membership problem can then be formulated as a 
linear program (LP) [21] 


maximize: A • p — C 

subject to: A • — C < 0 , for all (18) 

Ap-C^ 1. 


The variable of the LP are {A, C}, where A represents the coefficient of a Bell-Svetlichny inequality, detecting genuine 
multiparticle nonlocality, and C the corresponding local-nonsignalling bound. The LP optimizes the coefficients A to 
obtain the maximal value (at most C -I- 1) for the quantum probabilities, while keeping the local-nonsignaling bound 
C. As a consequence, the vector p can be written as a convex combination of {v^} if and only if the optimal value of 
the LP is 0. 

The noise tolerance for jiJa) can then be computed by mixing it with white noise, i.e., \H 3 ){H 3 \ i-Ag={l — 
e)\H){H\-\-el/8, and compute for which values of e the LP (18) gives optimal value 0. Standard numerical techniques 
for LP give that up to e Ri 1/13 Ri 7.69% the probabilities for X and Z measurements cannot be explained by a hybrid 
local-nonsignalling model. 
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APPENDIX B: CORRELATIONS FOR THREE-UNIFORM HYPERGRAPH STATES 


B.l. Preliminary calculations 


Before starting with the actual calculations, we need to settle couple of identities and a look-up table, which we 
will refer to throughout the main proofs. 

The first and probably the most important identity is a commutation relation between multi-qubit phase gates and 
Pauli X matrices [9], 

a((g)X,) = ((g)Xfc)( n ^eU/}). (19) 

i&K i&K f€V{K) 


Here, V{K) denotes the power set of the index set K. Note that the product of the CeVi/} ™ay include the term Cg, 
which is defined to be —1 and leads to a global sign. 

Furthermore, it turns out to be useful to recall some basic facts about binomial coefficients, as the appear frequently 
in the following calculations. 

Lemma 1. The following equalities hold: 


Re 


(i + i)' 


s : 

fc=0,4,... ^ 


n 

fc + 2 


Im 


4+*)"]= E 

k=0A,. 


n 

k + 1 


n 

k + 3 


Proof. Here we derive (20) and (21) together: 




n\ / ^ 

J+*U+i 


k 2/ i/c -|- 3 


k=o ^ '' fc=0,4, 

It is easy to spot that i?e[s] and Irnfs] indeed leads to the identities (20) and (21) respectively. 


( 20 ) 


( 21 ) 


( 22 ) 

□ 


The following look-up table represents the values of (20) and (21) for different n. These values can be derived from 
the basic properties of complex numbers: 


# 

n 

Re (1-fi)" 

Im (1 -I- z)" 

Re (1-1-z)” -\-Im (1-1-z)" 

Re (1-1-z)" —Im (1-1-z)” 

1. 

n = 0 mod 8 

-f2t 

0 

-k2t 

-f2t 

2. 

n = 1 mod 8 


+2^ 

+2^ 

0 

3. 

n = 2 mod 8 

0 

-k2t 

-k2^ 

-2t 

4. 

n = 3 mod 8 

-2^ 

+2^ 

0 

1 

to 

5. 

n = 4 mod 8 

-2t 

0 

-2t 

-2t 

6. 

n = 5 mod 8 

-2^ 

-2^ 

-2^ 

0 

7. 

n = 6 mod 8 

0 

-2t 

-2t 

-f2t 

8. 

n = 7 mod 8 

+2^ 

-2^ 

0 

-t2^ 


Table 0. Look-up table for the values of (20) and (21) . 


B.2. Correlations for X and Z measurements on fully-connected three-uniform hypergraph states 


Lemma 2. Consider an arbitrary number of qubits N and the three-uniform fully-connected hypergraph (HG) states. 
Then, if m is even with I < m < N the following equality holds 


{X ...XZ...Z) 


if m = 2 mod 4, 
— i if m = 0 mod 4. 


m 


(23) 
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Proof. We can write: 


K := {HG\X ...XZ ...Z\HG) = Ge) X ... X Z... ZI nC'e)|+)®^. (24) 

'^eeB ^ ^eeE ^ 

We can group all the controlled phase gates on the right hand side of the expression (24). Note that the operators Ce 
and Xi do not commute, but we can use the identity (19). While regrouping we count the multiplicity of each phase 
gate. If each phase gate appears even times, we get an identity as G^ = 1, if not, we keep these phase gates with the 
multiplicity one for the further calculations. 

For the purposes which will become apparent shortly, we denote the parties which measure in X direction by @ and 
ones in Z direction by A, in a way that, for example, if an arbitrary phase gate acts on XXXZZ^ it is represented 
as ® ® ® A A. Without loss of generality, we fix one phase gate Ce and consider all the possible scenarios of ® and 
A it can be acting on. Since we work on three-uniform HG states, every phase gate acts on bashes of different three 
party systems. These parties can be either of type ® or A and we have to consider all possible scenarios. Since we 
are working with symmetric states, we can sort the parties such that we have m ®’s followed by (TV — m) A’s: 



N—m 


and then we can enumerate all the scenarios of one phase gate acting on (25): 


1 . GeZZZ corresponds to A A A 3. GeXXZ corresponds to ® ® A 

2. GeXZZ corresponds to ® A A 4. GgXXX corresponds to @ ® ® 

We consider each case separately: 

1. GeZZZ = ZZZGe as Cg and Z commute. Cg moves on the right side with the multiplicity one. To save us 
writing in the future, we will denote the multiplicity of the phase gate moving on the right side by ffe. In this 
particular case it is # A A A = 1. However, on the right side of the equation (24) we have a product of all three-party 
phase gates. Therefore, we get Gg with the multiplicity of two and Gg = 1. Note that, as we have chosen an arbitrary 
three-qubit phase gate, the same result holds for every such phase gate. So, all three-qubit phase gates coming from 
the case 1, cancel out. We will see that with the same reasoning all three qubit phase gates cancel out (give an identity). 

2. For GgXZZ, we use the identity (19): 


GgXZZ = XCgG{g\x}ZZ. (26) 

The three-qubit phase gate Gg, from (26), appears with the multiplicity one {ff ® AA = 1) and like in the case 1, 
it gives an identity when being multiplied by the same three-qubit phase gate at the right side of the expression in 
(24). It is more tricky to calculate the multiplicity of G{g\x} (# A A as the ® part (or equivalently, X part) is 
removed from the set of vertices e.). For this we need to fix A A and count all the scenarios when an arbitrary Cg 
is reduced to AA. As we are working with the symmetric case, such scenario repeats (’(“) = m times, where m is 
the number of parties measuring in X direction.We shortly denote this as ^ A A = (™) = m. So, as m is an even 
number, {C{g\x})'^ = (Gaa)™ = 1 - 

Note that the gate Caa can only be generated from in case 2. 


3. For GgXXZ, we use the identity (19): 


GgXXZ = XXGgC{g\xx} [ n C'{e\x} 

VA 


Z. 


(27) 


The three-qubit phase gate Gg, from (26), appears with the multiplicity one (# ® ®A = 1); therefore, like in the two 
previous cases, it cancels out on the right side of the expression. A multiplicity of G[g\xx} is calculated by fixing a 



concrete A and counting all possible appearance of arbitrary ®@. As the number of parties measuring in direction is 
X is TO, this means that it is all combination of two parties with X measurements out of total to parties. So, 


fm\ to(to-I) 

-2- 

The last one from this case is the multiplicity of C'{e\jc} or ^ ® A. Here we fix one qubit from to {X direction) and 
one from N — m {Z direction) and count the number of such occurrences, when the third qubit is an arbitrary one of 
the type ®, which is exactly Therefore, 

/ ^ \ 

#®A= ]=m — 1, which is odd, C® a remains. (29) 


even, if to = 0 mod 4 
odd, if TO = 2 mod 4 


Ca cancels out. 
Ca remains. 


(28) 


4. For CeXXX, we use the identity (19): 


CeXXX = xxxc, 


[n 

\/x 


C. 


{e\X} 




C. 


{}■ 


(30) 


Ce occurs once and it gives an identity with the other one from the right side like in the previous cases. The multiplicity 
of C{e\x} is # ® @. Here we fix two parties in X direction and count the occurrence of this scenario by altering the 
third party, from the remaining to — 2, in X direction. Therefore, 

/ jYi, — 2 '\ 

#®®=( )=TO — 2, which is even, => C@@ cancels out. (31) 


Similarly for C{e\x:x}i we fix one party in X direction and count all possibilities of choosing two parties out of 
remaining to — 1. Therefore, 


/to — 1\ (to — 1)(to — 2) 

V 2 2 


odd, if TO = 0 mod 4 C@ remains, 
even if, to = 2 mod 4 C@ cancels out. 


(32) 


At last, we consider C{}. This gate determines the global sign of the expectation value and it appears only when Ce 
acts on systems which are all measured in X direction. Therefore, 


#{} = 


m{m — 1)(to — 2) 

a”3 


which is even => the global sign is positive . 


(33) 


To go on, we need to consider two cases #1: to = 0 mod 4 and #2: to = 2 mod 4 and calculate the expectation 
value separately for both: 

Case When to = 0 mod 4, we write out all remaining phase gates and continue the derivation from equation 
(24): 

(K) ■.= {+\®^X ...XZ...Z J]; C'@aC'@|+)®'^. (34) 

V®,VA 


Using the fact that X is an eigenstate of (+|, we can get rid of all Xs and then we can write Ca instead of Z-. 


^ 11...11 .. 11...11 

(x)=(+i®^ n ^^®ac'®cai+)®^=^ ^ (*i n 1^') 

v2 i=00...00 V®,VA v2 j=o 


V®.VA 


i=oo...oo 


1 

''W 


(00...00I Y[ C@aC@Ca|00...00) + --- + (11...11| Y[ C@aC@Ca|11...11) 

V®,VA V®,VA 


(35) 



C(g,ACig)CA 

-V®.VA 
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In (35), to get line two from the line one, note that Hv© VA C'@aC'@C'a is a diagonal matrix. 

To evaluate the trace of the given diagonal matrix, we need to find the difference between the number of +1 and 
— 1 on the diagonal. We write every row in the computational basis by enumerating it with the binary notation. For 
each row, we denote by a the number of I’s in binary notation appearing in the first m columns and by /3, the same 
on the rest. For example, for iV = 7 , and m = 4, the basis element jllOlllO) leads to a=3 and /? = 2. Considering 
the phase gates in the equation (35), the expression (—1)'* defines whether in the given row the diagonal element is 
+ 1 or —1, where : 

In s, (“) (^) denotes how many C^a acts on the row. Also, (“) determines the number of C® and (^), number of 
Ca- Every time when the phase gate acts, it changes the sign of the diagonal element on the row. Therefore, we 
need to determine the number s: 


To see whether s is even or odd, we have to consider the following cases exhaustively: 

1. a is even & /3 is even (—1)'* = +I 


2. a is even & /3 is odd (—1)'* = — I 

3. a is odd & /3 is even (—1)'’ = —1 

4. a is odd & /3 is odd (—1)'* = —1 

From the cases 3 and 4, one can directly calculate the trace: 

N- 


These two cases sum up to zero. 

These two cases contribute with the negative sign. 


{K) = 


I 


- E 

a=l,3,. 


E 

B=0 


N — m 




2N 


So, we get that if m is divisible by 4, 


(X_^Z...Z) = --. 


(37) 


(38) 


Case #2: We use the identical approach: when m = 2 mod 4, we write out all remaining phase gates and continue 
the derivation from equation (24): 


{K) = {+\®^X...XZ...Z CgACAlT) 

V®,VA 


0N 


(39) 


Again we use the fact that X is an eigenstate of (+| and Z = Ca- As in (40), there is already one (Ca), they cancel. 
Therefore, we are left with: 


(K) 


(+1®^ n = 

V®.VA 


C®A 


V@,VA 


(40) 


We need to define the sign of the diagonal element by (—1)'*, where 


s = 


a\ f fd 


1/ VI 


= ap. 


1. a is even (—1)® = +I 


This case contributes with the positive sign in the trace 
These two give zero contribution together. 


2. a is odd & /3 is even (—1)® = +I 

3. a is odd & /3 is odd (—1)® = —I 
As the case 2 and 3 add up to zero, we only consider the case 1: 




2N 


Q:=0,2,. 


N—m 

E 

/3=0 


N — m 
P 


(2m—l(2N—m 2 

^ “ 2' 


(41) 


(42) 
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So, we get that if m is NOT divisible by 4, 

m 


( 43 ) 


This completes the proof. 


□ 


B.3. Correlations for X measurements on fully-connected three-uniform hypergraph states 
Lemma 3. If every party makes a measurement in X direction, then the expectations value is 


{ XX XX) = 

N 


0 if N = 0 mod 4, 
1 if N = 2 mod 4. 


Proof. In this proof we employ the notation introduced in details in the proof of the lemma 2. 


{K) := {H^\XX ... XX\H^) = (+|®^ [ Ce 


XX ... XX 


e£E 




(44) 


(45) 


e£E 


We use the identity (19), to regroup the phase gates on the right hand side of the expression (45). Therefore, we 
count the multiplicity of the remaining phase gates: 

#® ® ® Each Ce, where |e| = 3, occurs once and cancels with the one on right hand side. 

#®® = (^ 7 ^) cancels 

^0 _ — {n-i){n- 2 ) .g I odd, if X = 0 mod 4 C® remains. 

^ ' 2 I even, if X = 2 mod 4 C@ cancels. 

~ (^) “ jg g.^gj^ ^ global sign GS is positive. 

Therefore, we need to consider two cases to continue the derivation of the expression (45): 

Case # 1 : If X = 0 mod 4, then 


(X) = (+|®^XX...XX Y[ = (+1®^ Y[ = 0. 

®&E ®&E 


Case ff2: If X = 2 mod 4, then 


(X) = (+|®^XX...XX|+)®^ = (+1+) 


N 


= 1 . 


(46) 

(47) 

□ 


B.4. Correlations for Z measurements on fully-connected three-uniform hypergraph states 

Lemma 4. If every party makes a measurement in Z direction, the expectation value is zero. Therefore, we need to 
show that 


(X) := {ZZ...ZZ) = 0 


(48) 


Proof. 


(X)=(+|®^[n^e 

eGB 


zz ...zz 


[lie. 

eeE 


| 0 )®Ar ^ , _ zZ\+)‘^^ = ±^Tr{Z ...ZZ)=0. 


(49) 

□ 
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APPENDIX C: PROOF OF OBSERVATION 5 

Observation 5. If we fix in the Bell operator Bn in Eq. (15) the measurements to be A = Z and B = X, then the N- 
qubit fully-connected three-uniform hypergraph state violates the classical bound, by an amount that grows exponentially 
with number of qubits, namely 

{Bn)c ^ for local HV models, and 

{Bn)o ^ \ for the hypergraph state. (50) 

V 2 

The classical bound can be computed using the bound for Mermin inequality, i.e., {B^)c < . If V is odd, our 

Bell inequality is exactly the Mermin inequality. If N is even, {Bn)c = ’ ^n-i + ^ ' ^n-i)c < 2(S)Jf_;^)c, where B 

denotes the inequality where A and B are exchanged, and the claim follows. The quantum value is computed directly. 
In particular, for odd N, we have {Bn)q = Y.k even If iV = 0 mod 4, then {Bn)q = Y.k~eLn {T)\ “ 

(A ... A) = 2^-2 - 1/2 and for A = 2 mod 4, we have {Bn)q = Ek eLn (fe) 5 “ = 2^”" + 1/2. 

Notice that, since Eq. (15) in the main text is a full-correlation Bell inequality, it can be written as a sum of 
probabilities minus a constant equal to the number of terms. With the proper normalization, such probabilities exactly 
correspond to the success probability of computing a Boolean function (cf. conclusive discussion and Ref. [29]). 


APPENDIX D: CORRELATIONS FOR FOUR-UNIFORM HYPERGRAPH STATES 

Lemma 5. The following statements hold for N-qubit, four-uniform hypergraph states: 

1. For the case N = 8k — 2 or N = 8k — 1, or 8k, we have: 

(*) 


{X ...XZ ...Z) = 


+ j/(m - 1) = 0 mod 4, 

if (m-l) = 2 modA. 


(ii) For N = 8k — 1, we have: 


( AA...AA ) = -1. 

N 


For N = Sk — 2 or N = Sk, these correlations will not he needed. 


(51) 


(52) 


2. For N = Ak -\-l, we have: 


{X ...XZ ...Z) = 


+ 2 rw 2 i : «/(m - 1) = 0 mod 4, 
- 2 rW 2 i : if {m - 1) = 2 mod A, 
if m = N. 


2l 

' 2 ^ 

1 

2L^/2J 5 


(53) 


3. For N = 8k -\- 2, or 8k -\- A, we have for even m: 

(i) 


{X...XZ...Z) = \*(SZ, ./(]V-™) = 0™<i4, 

'—''—^ // (-^ — m) =2 mod A. 


(54) 


(ii) 


{XX... XX) = 


2--^-hi 
2 ^ ■ 


(55) 
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4- For N = 8k + S, { X . ^^. Z \) for even m gives the same exact result as the part 3, so we have: 

m N—m—l 


1 -^m 7 ^ */(M-m) = 2 mod 4, 


( 56 ) 


where M = N — 1. 

Proof. Each part of the theorem needs a separate consideration, note that the notation and machinery that is 
employed in the proof is based on the proof of the Lemma 2. Therefore, we advise the reader to become familiar 
with that one first. 

Part 1. We consider the cases when N = 8k — 2, N = 8k — 1 or 8k and odd m together. 

We prove {i) first: 


(Gi) = {H^\X_^Z ... Z\H^) = (+1®"^ (^l[Ce)xX...XZ...z(^l[Ce) 1+)®"^. (57) 

m eeE e^E 

We need to use the identity (19) to regroup all the phase gates on the right hand side of the expression (57). If each 
phase gate occurs even number of times, they give an identity, otherwise, they are used in the further calculations. 
We consider each case separately in the following table: 


1. #A A A 

2. #AA 
#® A A 


is odd => Gaaa remains. 

=> Gaa cancels. 
Gaa remains. 

('"("^) is even G@aa cancels. 


/m\ m(m—1) 

V2/ “ 2 


IS 


even, if (m — 1) = 0 mod 4 
odd, if (m — 1) = 2 mod 4 


3. #A 

#® @ A 
#®A 


4. #® ® @ 

#®® 

#{} 


(T) = 

(^-2) 

(^- 3 ) 


n(?7i—l)(?7i—2) 

2^3 


IS 


even, if (m — 
odd, if (m — 


is odd => G@@a remains. 


(m—l)(m—2) 


IS 


even, if (m — 
odd, if (m — 


1) = 0 mod 4 
1) = 2 mod 4 

- 1) = 0 mod 4 
1) = 2 mod 4 


=> Ga cancels. 

=> Ga remains. 

G@a cancels. 
=> G@a remains. 


is even => G®®® cancels. 

_ (m- 2 )(m- 3 ) J odd, if (m - 1) = 0 mod 4 ^ G®® 

IS \ 

I even, if (m — 1) = 2 mod 4 G®® 

_ (m-i)(m- 2 )(m-^ even => G® cancels. 
(m-i)(m-2Xm-3)(m-4) global sign (GL) and will 

Table 1. Counting phase gates for a four-uniform case. 


remains. 

cancels. 

be discussed separately, 
m is odd. 


Remark: All four-qubit phase gates move with the multiplicity one to the right side and therefore, cancel out with 
the same phase gate on the right. The detailed reasoning was discussed in proof of the Lemma (2) in the Appendix 
A. So, we skipped such scenarios in the Table 1. 

Now we consider two cases of (m— 1) separately and for each case we fix the global sign (GL) defined in the Table 1. 
Case #1: (m — 1) = 0 mod 4: 

{Gi)=±{+\®^XX...XZ...Z n GAAAG®®AG®®|+)®'^ = ±^rr( n GaaaG®®aG®®Ga). (58) 

V@,VA V@,VA 

Remark: We write the ’±’ sign as we have not fixed the global sign yet. 
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To evaluate the trace of the given diagonal matrix, we need to find the difference between the number of +l’s 
and — I’s on the diagonal. We write every row in the computational basis by enumerating it with the binary 
notation. Due to the symmetry of the problem, we assign the first m columns to X measurement (®) and the rest 
to, Z{A). For each row, we denote by a the number of I’s in binary notation appearing in the first m column 
and by /3, the same on the rest. This notation is also adopted. See the proof of Lemma 2 for more detailed explanation. 

Considering the phase gates in (58), the expression (—1)® defines whether in the given row the diagonal element is 
+1 or —1, where : 

--(:.)-(;) (0 - (3 ^ 0 ^ . 1 ).. 

The sign of the diagonal element is determined at follows: 

1. a is even & /3 is even: if a = 0 mod 4 (—1)'* = +1 

if a = 2 mod 4 => (—1)'’ = —1 

2 . a is odd & /? is even: if (a — 1) = 0 mod 4 => (—1)® = +1 

if 0 — 1) = 2 mod 4 => 01)® = —1 

3. (Any a) & ,5 is odd: if (/3 — 1) = 0 mod 4 (—1)® = —1 

if 0 — 1) = 2 mod 4 (—1)® = +1 

Having established the ±1 values for each row, we can sum them up to find the trace in (58). Here we use the 
identities (20) and (21) and afterwards the look-up Table 0 to insert the numerical values where necessary: 


(Gi) = ± ^ 


+ E 

/3=1,5,. 


N—m 

E 

/3=0,2,4... 


N — m 

p 


E 

^ a. — 0,4:,. 


m \ f m 

a-\-l) Va -h 2 


m 

G; -|- 3 


N — m\ (m\ f N — m\ 

E(„) + ( ,3 + 2 )E 


= ± 


1 


^ 2^ 


Re 


Re 


(i+*r 

(i+*r 


■ Im 


■ Im 


(i+*r 

(i+*r 




N—m r 

1 + 2™ 

13=1,5,... 


N — m\ f N — m 
P + 2 


^N-m-l _ 2m 


(1+0 


N-r. 


= +—at-G. 

2N 


We have to consider N = 8k — 1 and N = 8k or N = 8k — 2 separately to continue the derivation of (60): 
1. For N = 8k — 1, using the values from the Table 0: 




2^2^-™-! -2™/m 


(1 + 0 


N—m 


^2^2^“™“^ -I- 2™2’^V^ + 


2N 


2W/2J-Lm/2J ■ 


(60) 


(61) 


2. For N = 8k orN = 8fc — 2 , using the values from the Table 0: 


(Gi) = ± 


Therefore, 


2N 


2^2^-™-^ _ 


(1 + 0 


N—m 


= ±- 


2^2^-™-! , 2m2 ^ T ^ 


2N 


= ± 


2 W/2j-m _|_ Y 
2W/2J-Lm/2J ■ 


(Gl) - ±2[Ar/2j-Lm/2j ■ 


(62) 


(63) 


Concerning the sign in (63), it is affected by the product of two components: one from the case 4 from Table 1: 
GL and the other by E from (60). If (to — 1) = 0 mod 8, the equation E has a positive sign and also GL = -|-1. 
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And if (m — 1) = 4 mod 8, E has a negative sign and GL = —1. Therefore, in both cases or equivalently, for 
(to — 1) = 0 mod 4, 


{Gl) - + 2[N/2\-[m/2i ■ 


(64) 


Case #2: (to — 1) = 2 mod 4: 


{Gl) =±{+\^^XX...XZ...Z n CaaaC'aaC'®®aC'@aCa|+)®'^ 

V®,VA 

= C'aaaC'aaC'@@aC'@a)- 

V®.VA 


(65) 


In this case, we apply the same technique to determine the sign: 


s := 



The sign of s, is determined at follows: 


/a\ //3\ ^ /3(/3-l)(/3-2) /3(/3 - 1) a{a - 1) 

2-3 2 2 


1. /3 is even & any a if /3 = 0 mod 4 => (—1)^ = +1 

if /3 = 2 mod 4 => (—1)^ = —1 

2 . /3 is odd & a is even: if a = 0 mod 4 => (—1)® = +1 

if a = 2 mod 4 => (—1)® = —1 

3. (3 is odd & a is odd: if (a — 1) = 0 mod 4 => (—1)® = —1 

if (a — 1) = 2 mod 4 (—1)® = +1 


( 66 ) 


(Gi)=± 


1 

w 

E 

= 1,3,. 
1 

w 



= ± 


2™ E 

/3=0,4,. 


N — m\ f N — m 
i + 2 


■\N—m — 1 


Re 


2^Re 


(1 + 0 


N-r, 


+ 2 


N—m—l 


Re 


(1 + 0^ 


— Jrn 


TO 

0 + 3 

(1 + 0 ’ 

(1 + 0 "’ 


— Im 


(1 + 0 ^ 


- 2N^- 


We have to consider iV = 8fc — 1 and N = 8k or N = 8k — 2 separately to continue the derivation of (60): 
1. For N = 8k — 1, using the values from the Table 0: 

1 


(Gl) = ± 


2N 


2™i?e 


(1 + 0 


N-r 




^ 2 '^ 2 ^^ + 2 ^ 2 ^-'^-^ ^2L^/2J-”^ + 1 


2N 


2W/2J-Lm/2J ■ 


2. For N = 8k orN = 8k — 2 , using the values from the Table 0: 


(Gi)=±^ 


Therefore, 


2™i?e 


(1 + 0 


N—m 






2N 


(Gl) = ± 


2LAf/2j-m _|_ 2 ^ 


(67) 


■2LW2J-Lm/21 ■ 


( 68 ) 


(69) 


2LW/2J-Lm/2J ■ 


(70) 
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Concerning the sign in (70), it is affected by the product of two components: one from the case 4 from Table 1: GL 
and the other by E in (67). 

Hence, if (m — 3) = 0 mod 8, E has a negative sign and GL = +1. And if (m — 3) = 4 mod 8, E has a positive sign 
and GL = —1. Therefore, in both cases or equivalently, for (m — 1) = 2 mod 4, 

2\_N/2\-m ^ 

^ “ 2VN/2\-ym/2\ ■ 

This completes the proof of part 1 (i). 

Part 1 {ii) For iV = 8fc — 1, show that 

(G 2 ) = ( AA...AA ) = -1. (72) 

N 

Here as well we use the identity (19) to count the multiplicity of remaining phase gates. Since all the measurements 
are in X direction, we need to make a new table with the same notations as in the previous case: 


1 . 


every gate occurs only once => every Gg cancels with the Gg on the right hand side. 


2 . 


(8.-4) 


is even cancels. 

3. # ® ® (*^ 2 "^) ^ cancels. 

4. (*^ 3 "^) is even => G® cancels. 

5. #{} is odd => we get a global negative sign. 

Table 2. Counting phase gates for a four-uniform HG when each system is measured in X direction. 


Therefore, 


(G2) = -^Tr(l)=-l. 


(73) 


This finishes the proof of part 1. 

Part 2: We show that, for = 4fc -|- 1: 


- 2 \ m/ 2 \ 1 if (to — 1) = 0 mod 4, 
(G 3 ) = (^._XZ ■■■Z)={ if (to - 1) = 2 mod 4, 

2 TW> \im = N. 


(74) 


Since the number of systems measured in the X direction is the same in this part as it was in part 1, we can 
use the results demonstrated in the Table 1. Therefore, we use equation (60) when m — 1 = 0 mod 4 and (67), for 
TO — 1 = 2 mod 4. 


Case #1: (to — 1) = 0 mod 4: 

1 


(G3) = ± 


2N 


Re 


(i+*r 


Im 


(i+*r 


2^-'"-! - 2'"/to 


(1 + 0 


N—m 


= ±^E. 

2N 


As = 4A:-|-1, we have that A^—m = 4fc-|-l —to = 4fc—( to— 1 ), which is divisible by 4. Therefore /to (1-l-i)'^ 
and equation (75) reduces to: 

2^N—m — l 


(G 3 ) = ±- 


2N 


Re 


(i+*r 


Im 


(i+*r 


= ± 


1 

2r™/2i 


(75) 


= 0 . 


(76) 
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We need to fix the global sign GL from the Table 1. For this, we consider two cases. First, if (m — 1) = 0 mod 8, 
then GL = + and so is the sign E in equation (75): 




(77) 


Second, if (m — 1) = 4 mod 8, then GL = — and so is the sign ofif in equation (75): 


(Ga) 


1 


(78) 


Case # 2: For (m — 1) = 2 mod 4: 


(Ga) = 


2^Re 


(1 + 0 


iV-ri 


'iN—m — 1 


Re 


(i+*r 


— Im 


(i+*r 


(79) 


As N = 4fc + 1, we have N — m = Ak + l — m = 4:k— (m — 1), which is not divisible by 4 but is an even number. 


Therefore, Re 


( 1 - 


\N—m 


= 0. So, the equation (79) reduces to: 


(Ga) = ±- 


-yN—m—1 ' 


2N 


Re 


(i+*r 


— Im 


(i+*r 


■yN—m—l 


= ±- 


2N 


-E. 


(80) 


We need to fix the global sign GL from the Table 1. For this, we consider two cases. First, if (m — 3) = 0 mod 8, 
then the global sign is positive but the sign of E in (80) is negative. Therefore, 


(Ga) = 


1 


2rm/2l ■ 

Second, if (m — 3) = 4 mod 8, then the global sign is negative but the sign of E in (80) is positive. Therefore, 

1 


(Ga) = - 


2rm/2l 


(81) 


(82) 


Case #3: m = N resembles part 1 (ii). The only difference comes in with the number of qubits we are currently 
working with: 


1 

2 . 

3. 

4. 

5. 


# ® ® ® @ e every gate occurs only once every Ce cancels with the Ge on the right hand side. 

# ® @® is even => C®®® cancels. 

# ® @ (^^ 2 "^) ^ G@® remains. 

#@ is even => C® cancels. 

#{} the global sign depends on fc, as = ("^^+i)4fc(4fc-i)(4fc-2) 

Table 3. Counting phase gates for a four-uniform HG when each system is measured in X direction. 


Back to the expectation value. 


(G3) = ±^Tr 


\ G®® 

V® 


(83) 


So, we have to count the difference between the amount of +l’s and — I’s on the diagonal. As we use exactly the 
same techniques before, we will skip the detailed explanation. The sign on the diagonal is: 


(- 1 )( 5 )=(- 1 )^ 


( 84 ) 
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and it is straightforward to evaluate it for each value of a. 


(Ga) = ± 


2N 


N 

E 

q;=0,4,... 


N' 


+ 


N 

a + 1 


N 

a + 2 


N 

cr “t" 3 


= ± 


2N 


Re 


( 1 +*) 


N 


-\-Im 


( 1 +*) 


N 


= ±^E. (85) 


Keeping in mind that = 4fc + 1, the global sign from the Table 3 is positive for even k and negative for odd. The 
the sign of E in (85) is positive if k is even and negative, otherwise. Therefore, 

^^ 3 ) = 2[Ar/2j ■ 

This completes the proof of part 2. 

Part 3: We start with (ii). We show that for TV = 8 A: + 2, or 8 fc + 4: 

(ii) 

2^“i + 1 

{Gi) = {XX...XX) = - (87) 

2 2 

Although the result seems identical, unfortunately, each case needs a separate treatment. The technique is similar 
to the previous proofs, though. We just mind the number of qubits we are working with: 

For = 8 A: + 2 we find the remaining phase gates as follows: 


1 # (i) @ ® @ each gate only once; thus, every Ce cancels with the Ce on the right hand side. 

2 . # ® @® is odd => C®®® remains. 

3. # ® @ (® 2 ^) is even => C®® cancels. 

4. #@ is even => C® cancels. 

5. #{} (*^ 4 *’^) is even times => we get a global positive sign. 

Table 4. Counting phase gates for a four-uniform HG when each system is measured in X direction. 


Therefore, 


(G4) — ^T’J'I^J^G®®® 

V® 


( 88 ) 


We use (—1)® to define the sign of the diagonal element and s = ( 3 ). So, after considering all possible values of a, 
it is directly obtained that 


(G4) 2Af^^(G®®®) 2N 


1 


N 


E 

q;=0,2,... 


E 

0=1,5,... 


r/N 
a 


N 

a + 2 


(89) 


■)N-1 


/ to [(1 


\N-\ 


2 ^- 1 -kl 

2 f ■ 


For N = 8 fc -I- 4 we find the remaining phase gates as follows: 
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1 - # ® @ ® @ each gate occurs only once => every Ce cancels with the Ce on the right hand side. 

2. # ® ®® is odd ^ C@@@ remains. 

3. # ® @ (*^ 2 *"^) ^ remains. 

4. #@ is odd = C@ remains. 

5. #{} is odd we get a global negative sign, GL = —1. 


Therefore, 


Table 5. Counting phase gates for a four-uniform case in all X direction. 


(G4> = -2^rr[n 




(90) 


We use (—1)® to define the sign of the diagonal element and s = ( 3 ) + ( 2 ) + (i)- So, after considering all possible 
values of a, it is directly obtained that 


{G 4 ) — — 


N 


E 

Q!=0,4,... 


N 

a + 2ji 


1 

'W 




N 

- E 

0=1,3,... 

2^-1 + 2^/2 


(91) 


2 --^ + 1 
2 f ■ 


This finishes the proof of part {i). 
(ii) We need to show that 


{G,) = {X...XZ...Z) = \ +3^ it (W - ,n) = 0 4, 

'—''—^ \ — ^ 2^/2 if (-^ ~ ™) = 2 mod 4. 


(92) 


Note that in this case m is an even number. Therefore, we have to derive again from the scratch how phase gates can 
be moved to the right hand side of the expression and for this we use the identity (19). 


1. #A A A 

2. #AA 
#® A A 

3. #A 

#® @ A 

#@A 

4. #@ ® @ 

#®@ 

#{} 


(*(*) is even => Caaa cancels. 


( 2 ) = 


m(m—1) 


IS 


even, if m = 0 mod 4 
odd, if TO = 2 mod 4 


■ Caa cancels. 
Caa remains. 


^) is odd => C(gAA remains. 


(T) = 

(m—2 


_ m(m—l)(m— 2 ) 


2-3 


is even => Ca cancels. 


(™^ ^) is even => C^^a cancels. 


— 1 \ (m—l)(m— 2 ) 


('" 2 -^) = 


IS 


odd, if TO = 0 mod 4 
even, if to = 2 mod 4 


C(gA remains. 
■ C(gA cancels. 


("*-3) is odd ^ Cg 

{m — 2\ _ (m—2)(m—3) 

V 2 / ~ 2 


IS 


y®®® remains. 

odd, if TO = 0 mod 4 C(g@ remains, 
even, if to = 2 mod 4 => C®® cancels. 

odd, if TO = 0 mod 4 => C® remains, 
even, if to = 2 mod 4 C® cancels. 

(J2) = affects the global sign (GL) and will be discussed separately. 

Table 6. Counting phase gates for a four-uniform HG state, for even to. 


m — l\ (m—l)(m—2)(m—3) 


(A’’) = 


2-3 


IS 
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Remark: Similarly to previous proofs the four-qubit phase gates cancel out. Therefore, we directly skip the discussion 
about them. 

We need to consider two cases, when m = 0 mod 4 and m = 2 mod 4 for each TV = 8 fc -I- 2 and 8 fc -|- 4 separately: 
Case ^ 1: If m = 0 mod 4: 


{Gi)=±{+\®^XX...XZ...Z [] 

V@.VA 


= 

2N 


C@aaC'@aC®@®C'@@C’@C'a 

V@,VA 


(93) 


We use (—1)® to define the sign of the diagonal element and s = (“) ( 2 ) + (“) (^) -I- ( 3 ) -I- ( 2 ) + (“) + (f)- If s is 
even, the value on the diagonal is -1-1 and —1, otherwise. We consider all possible values of a and /3: 


1. a is even & /3 is even if a = 0 mod 4 (—1)® = +1 

if a = 2 mod 4 (—1)* = —1 


2. a is even & /3 is odd: if a = 0 mod 4 => (—1)® = —1 

if a = 2 mod 4 => (—1)® = -1-1 

From here one can easily spot that for even a, there is equal number of -|-1 and —1 on the diagonal. So, they do 
not contribute in the calculations. We now consider the odd a: 


3. a is odd & /? is even if /5 = 0 mod 4 => (—1)® = —1 

if /3 = 2 mod 4 => (—1)® = -1-1 


4. a is odd & /? is odd: if (/5 — 1) = 0 mod 4 (—I)® = —1 

if (/5 — 1) = 2 mod 4 (—1)® = +1 


We now continue calculation of the trace from (136): 


1 


(*^ 4 ) = 


= ±- 


2iv 

2 


a=l,3,5. 
m —1 

— Re 


N—m 

E 

9=0,4,.. 


2N 


(1 + *) 


N-r 



/TV — m\ /TV — rri\ /TV — m\ 
V/3 + 1 y ^ V/3 + 2 j ^ U + 3 I 


Im 


(1 + 0 


N-r. 


= ±- 


2N 



(94) 


We have to take care of the sign which appears from the product of the sign of the sum of real and imaginary 
part in (94) and global sign {GL), which we defined while deriving the remaining phase gates. If m is divisible by 8 , 
GL = -1-1 and since we are in TV = 8 /c + 2 case, (TV — m = 8k + 2 — m) — 2 mod 8 and therefore: 


(G 4 ) 



2 '//i /2 —1 
2^/2 ■ 


(95) 


If m is not divisible by 8 , the global sign GL = — 1, and the sum of real and imaginary part also contribute with a 
negative sign. Thus, 


(G 4 ) 


2N 



2^ml2—l 

2 W 2 


(96) 


Since the TV = 8 /c -I- 4 case is identical, we only have to mind the sign of the sum of the real and imaginary part. 
Here as well we consider two cases: if m is divisible by 8 , then the global sign GL = -|-I, and the sign of the sum of 
real and imaginary part is ” — ”. Therefore, 


(G4) 


2N 



2m/2—1 

+ 


2N/2 


(97) 
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And if m is not divisible by 8 , GL = — 1, and the sign of real and imaginary part is ” + ”. Therefore, 




2 ^/ 2 —! 

2^/2 


(98) 


Case # 2: If m = 2 mod 4: 

{Gi) = ±{+\®^XX...XZ...Z n C@aaC'@@@C'aa|+)®'^ 

V®,VA 

1 (99) 

— C'@aaC'®@@C'aaCa)- 

V®,VA 


We use (—1)'* to define the sign of the diagonal element and s = (“) ( 2 ) + ( 3 ) + ( 2 ) + (i)- If s is even, the value on 
the diagonal is +1 and —1, otherwise. We consider all possible values of a and /3: 

Considering the terms from odd a: 

1 . a is odd & /? is even if (a — 1) = 0 mod 4 => (—1)® = +1 

if (a — 1) = 2 mod 4 => (—1)® = —1 

2. a is odd & /? is odd: if (a — 1) = 0 mod 4 => (—1)® = —1 

if (a — 1) = 2 mod 4 => (—1)® = +1 
It is easy to see that these cases adds up to 0. 

3. a is even & /3 is even if /3 = 0 mod 4 (—1)® = +1 

if/3 = 2 mod 4 ^ (-1)® = -1 

4. a is even & /3 is odd: if (/3 — 1) = 0 mod 4 => (—I)'* = — 1 

if (/3 - 1) = 2 mod 4 ^ (-1)® = +1 

Therefore, 


(G4) = ± 


1 

w 


E 


a=0,2,4. 


= ±2 


m —1 


N—m 

E 

3 = 0 , 4 ,.. 

Re 


N — m 


N — m 

/3 + 1 


N — m 
13 + 2 


N — m 
' + 3 


(1 + 0 


N-r, 


— Im 


(1 


\N-r. 


= ib- 


2^Ti/2 —1 
2^/2 • 


( 100 ) 


To fix the sign, we need to first consider = 8 A: + 2, and (to — 2) = 0 mod 8 . Then the global sign GL = +1 and 
type of IV — TO also yields a positive sign. But if (to — 2) = 4 mod 8 , global sign in negative and the N — m also 
yields the negative sign. So, 


(G 4 ) = 


2^ti/2—1 

2^/2 ■ 


N = 8k + A case is identical to = 8 A: + 2, therefore we will just state the result. For N = 8k + A 

2mj2—l 


(G4) = - 


2 N /2 


( 101 ) 


( 102 ) 


To sum up. 


(G4) = 


+ ^ 2 iv /2 if (A^ — to) = 0 mod 4, 


.1/2-1 


— ^^jv /2 if {N — to) = 2 mod A. 


(103) 


This finishes the proof of part 3. 
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Part 4: We show that for N = 8fc + 3, ( X . ^■ y even m gives the same exact result as the part 3 (i). 

m N—m—1 

We tackle the problem as follows: We make a measurement on one of the qubits in Z direction and depending on 
the measurement outcome, we obtain the new \H^ ) state, where M := fV — 1. Then we consider the expecta¬ 
tion values for the all possible measurement outcomes and From that we conclude the statement in the part 4. 


Initial HG state \H^) can be written in the following form as well: 


(104) 


where e' represents the gates containing the last qubit, N and e" represents the ones which do not contain qubit. 
And e = e' -|- e". Then if one makes a measurement in Z basis on the last qubit and obtains outcome -I-, 


= {0n\H^) = [(0|Arle/\Ar|0)Ar(0|Ar -I- (0|ArC'e'\Ar|l)Af(l|Af]C'e"|+)®^ 

e' ,e" 

= le'\Ar(0|ArC'e"|-|-)®^ = (OIat le'yArCe" |+)® ^ (10) AT -|- 11) AT) ^ 

e',e" e' ,e" 

= n ^e'\NCe"\+r^ = X{C,"\+r^■ 

e',e" e" 


(105) 


So, -b outcome after measuring in Z direction leaves us with which is precisely four uniform M-qubit HG 

state. Now, let us see what is the remaining state if one gets — as an outcome result: 


= {InW^) = [(l|7vle'\Ar|0)Ar(0|Ar + (l|ArGe/\Ar|l)Af(l|Af]C'e"l+)®^ 

e',e" 

“ n *^e'\Af(l|AfC'e" 1+)®^ = (l|w Ce'\NCe"\+)®^ {\0) N + |l)Ar) 
e',e" e',e" 

= n Ce>\NCe»\+)<^^ ■ 

e' ,e" 


(106) 


So, — outcome after measuring in Z direction leaves us with which is precisely a symmetric M-qubit HG 

state with all possible edges of cardinality four and three. We will call such HG state a three- and four-uniform HG 
state. 


Therefore, problem boils down to showing that, 

(i) If the measurement outcome is -b, we get the M = 8fc + 2 four-uniform HG state and the correlations are given in 
part 3. 

(m) If the measurement outcome is —, we get M = 8k + 2 three- and fouruniform HG state and the following holds: 


(Gs”) = {Hr- \^.^Z... Z\Hr:) = 


— if (A1 — to) = 0 mod 4, 

-b ^ 2 M /2 if — to) = 2 mod 4. 


(i) where M = 8k + 2 was already considered in part 3. 

(ii) . For \Hr-), 


(G-) = (+1®^ 


n Ge'Ce 

e' ,e"£E 


X ...xz ...z 


e' ,e"£E 


Before, we treated three- and four-uniform cases separately. Now, we just need to put them together. 


(107) 


(108) 


Case # 1: If m = 0 mod 4: 

Then from equations (35) and (136), we can directly write down that 


(GP = 


'-V®,VA 


(109) 










22 


We check the sign of each term on the diagonal by (—1)®, where s = (“) ( 2 ) + ( 3 ) + ( 2 ) + (f)- For this we need to 
consider each value of a and /? separately. 


1 . 


a is even & /3 is even 


a = 0 mod 4 => +1 
a = 2 mod 4 —1 


2 . 


a is even & /3 is odd 


a = 0 mod 4 => — 1 
a = 2 mod 4 => +1 


These two give zero contribution together. 


3. 

4. 


a is odd & /? is even 


a is odd & /? is odd 


/3 = 0 mod 4 +1 

/3 = 2 mod 4 => —1 

/3 — 1 = 0 mod 4 => — 1 
/3 — 1 = 2 mod 4 +1 


(GP = E 

a odd 


= ±- 


2 M 


M—m / -. r 

f M -m 

'-^= 0,4 

\M- 


M — m 

/3 + 1 


M — m 
f3 + 2 


M — m 
d + 3 


+ - Im{l + i) 


M-r. 


= ± 


2 ^ 


2 M /2 ■ 


( 110 ) 


If m = 0 mod 8 , real and imaginary part in (110) has a negative sign and the global sign coming from Table 6 , 
GL is positive. Note from equation (33) that three uniform gate moving does not introduce any global signs. And if 
m = 4 mod 8 , real and imaginary part in (110) has a positive sign and the global sign coming from Table 6 , GL is 
negative. Therefore, 

= ( 111 ) 

Case # 2: If m = 2 mod 4: 

Then from equations (40) and (121), we can directly write down that 


{G-,)=±^Tr 


C@aaC'@aC'@®®C'aa 

V®,VA 


( 112 ) 


We check the sign of each term on the diagonal by (—1)®, where s = (“) ( 2 ) + ( 3 ) + ( 2 ) + (“) (f) • For this we need 
to consider each value of a and /3 separately. 


1 . 

2 . 

3. 

4. 


a is even & ,5 is even 


= 0 mod 4^+1 
= 2 mod 4^-1 


a is even & /3 is odd 


/3 

/3 


1 = 0 mod 4 => +1 
1 = 2 mod 4 — 1 


a is odd & /? is even 
a is odd & /? is odd 


a — 1 = 0 mod 4 
a — 1 = 2 mod 4 => 
a — 1 = 0 mod 4 => 
a — 1 = 2 mod 4 => 


+ 1 

I => These two give zero contribution together. 

+1 


m = ± 2 M E 




M — m\ /M — m\ 

, /3 + 2 j ■ (v /3 + 3 )_ 


I ^ pm—1 

2 ^ 


Re{l + + Im{l + i) 


M—m 


= ± 


2 t-i 

2 M /2 ■ 


(113) 


If m — 2 = 0 mod 8 , real and imaginary part in (122) has a positive sign and the global sign coming from Table 
6 , GL is positive. Note from equation (33) that the three uniform gate moving does not introduce any global signs. 
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And if m — 2 = 4 mod 8, real and imaginary part in (122) has a negative sign and the global sign coming from Table 
6, GL is negative.Therefore, 


(G 5 -) 


+ 


2f-i 
2^/2 ■ 


Finally, we can put everything together. Since one can observe that {G^ ) 


-{Gt), 


|0)(0|(G+) - |1)(1KG5-) = |0)(0|(G5+) + |1)(1|(G5+) = 1(G+) = 

m N—m — 1 


This completes the proof of part 4 and entire lemma. 


(114) 


(115) 


□ 


APPENDIX E: BELL INEQUALITY VIOLATIONS FOR FULLY-CONNECTED FOUR-UNIFORM 

HYPERGRAPH STATES 


Lemma 6. An arbitrary four-uniform HG state violates the classieal bound by an amount that grows exponentially 
along with number of qubits. 

Proof. To show this, either in Eq. (15) or in the original Mermin operator, we fix A = Z and B = X. The 

choice of the Bell operator depends on the number of qubits: From Lemma (5) for a given N either the correlations 
for an even m or an odd m are given. If m is even, we choose Eq. (15) and (BfiJ), otherwise. 

From Lemma 5, it is clear that we need to consider separate cases. However, we choose the one which encounters 
the smallest growth in violation and this is the iV = 8fe + 3 case. Other cases just encounter different factors in the 
growth or are greater. For the Y = 8fc + 3, the strategy consists of measuring the Bell operator from Eq. (15) on 
M = N — 1 qubits. Then we have: 


{Bn)q > 



/ 1 \Af+2 


(116) 


Checking the ratio of the quantum and classical values, we have that 




{Bn)q 

{Bn)c 


N—^oo 4 


N-1 




■Ar+3 


1.20711^ 
2^2 + 2 ■ 


(117) 


Looking at the expectation values from Lemma 5, it is straightforward to see that in all other cases of N, correlations 
are stronger than in the N = 8fc + 3 case, so the quantum violation increases. □ 


APPENDIX F: BELL AND SEPARABILITY INEQUALITY VIOLATIONS FOR FULLY-CONNECTED 
FOUR-UNIFORM HYPERGRAPH STATES AFTER LOOSING ONE QUBIT 


Lemma 7. The following statement holds for Y = 8fc + 4 qubit, four-uniform complete hypergraph states: 


(Ge) = (Y_^_YZ_^1) 

m N—m — 1 



N—m+2 


N—m+2 


if m = 0 mod 4, 
if m = 2 mod 4. 


(118) 


Proof. The derivation of this result is very similar to the combinatorial calculations in the appendices D and E. Since 
m is even, we refer to the Table 6 to see what gates remain after regrouping hyperedges on the right hand side of the 
expression (118): 
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Case # 1: If m = 0 mod 4: 


(Gq) — ^ . ■ . Z 1 C@aaC'@aC'@a<c>C'@<c>C'@@@C'@0 C@ 1+)®’^ 

m N—m — l V@,A,^ 


= ±—Tr 
2 ^ 


C@AAC’®AC®A<0>C’®<C>C'@@@C'(g@C'@CA 

V@,A,^ 


( 119 ) 


Here @ again refers to X operator, A to Z and <C> to 1 and it is denoted as 7 . The strategy now is similar to the 

previous proofs: count the number of +l’s and -I’s on the diagonal and then their difference divided by 2^ gives the 

trace. 

We use (—1)® to define the sign of the diagonal element and s = (“) ( 2 ) + (“) (f) + (“) (f) (^) + (“) (^) + ( 3 ) + ( 2 ) + 

(?)+(?). It s is even, the value on the diagonal is +1 and —1, otherwise. We consider all possible values of a, /3, and 7 : 

a) If 7 is even (that is 7 = 0): 

1 . a is even & /3 is even if a = 0 mod 4 => (—1)® = +1 

if a = 2 mod 4 => (—1)® = —1 

2. a is even & /3 is odd: if a = 0 mod 4 => (—1)® = —1 

if a = 2 mod 4 => (—1)® = +1 

From here one can easily spot that for even a, there is equal number of +1 and —1 on the diagonal. So, they do 
not contribute in the calculations. We now consider the odd a: 

3. a is odd & /? is even if ^ = 0 mod 4 =7 (—1)® = —1 

if ^ = 2 mod 4 =7 (—1)® = +1 

4. a is odd & /3 is odd: if (/3 — 1) = 0 mod 4 =7 (—1)® = —1 

if (/3 — 1) = 2 mod 4 =7 (—1)® = +1 

b) If 7 is odd (that is 7 = 1): 


The cases 1 and 2 don’t change. Therefore, they sum up to 0. 

3. a is odd & /? is even if /3 = 0 mod 4 => (—1)® = +1 

if /3 = 2 mod 4 => (—1)® = —1 


4. Stays the same as in the previous case. 
Therefore the result is: 


(Ge) =±^( 


E 


E 

7 odd 


E 

OL odd 


odd 

r N—m—l 

E 

f=0,4... 


r N—m — l 

E - 

:=0.4... 


N — m — l 


N — m—l 
13 + 1 


+ 


N — m—l 
13 + 2 


N — m—l 
/3 + 3 


iV-m- 1 


N — m—l 
/9 + 1 


N — m—l 

/3 + 2 


iV — TO — 1 

/3 + 3 




1 

71 


N-m-\-2 


( 120 ) 


It is time to fix a sign. One needs to keep in mind that the sign of the Eq. (120) is negative: if to = 4 mod 8, the 
global sign from the Table 6 is negative, and /to[(1 + 1 )^“™“^] = —2 2 . Therefore, an overall sign in negative. 

If TO = 0 mod 8, global sign is positive and /to[(1 + = 2 2 . Therefore, an overall sign in negative. 
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Case # 2: If m = 2 mod 4: 

(Ge) = Y[ G@aaG@@@GaaG@a^Ga^|+)®'^ 

m A/"—m—1 V®,A,<0> 

=±^i’''( n G0AAG@@@GAAG(gA<>GAoGA)- 

V@,A,0 


( 121 ) 


We use (—1)® to define the sign of the diagonal element and s = (“) ( 2 ) + (“) (^) (1() + ( 3 ) + ( 2 ) + ( 3 ) (J) + ( 3 )- If 
s is even, the value on the diagonal is +1 and —1, otherwise. We consider all possible values of a, /3, and 7 : 


a) If 7 is even (that is 7 = 0): 

Considering the terms from even a: 

1. a is even & ,5 is even if /3 = 0 mod 4 => (—1)® = +1 

if /3 = 2 mod 4 => (—1)® = —1 


2. a is even & /3 is odd; if (/3 — 1) = 0 mod 4 => (—1)® = —1 

if (/3 - 1) = 2 mod 4 ^ (-1)® = +1 


Considering odd a: 

3. a is odd & /? is even if (a — 1) = 0 mod 4 => 

if (a — 1) = 2 mod 4 => 


4. a is odd & /3 is odd: if (a — 1) = 0 mod 4 => 

if (a — 1) = 2 mod 4 => 

It is easy to see that cases 3 and 4 adds up to 0. 


( 

( 

( 

( 


1 )® = +1 
1 )® = -1 

1 )® = -1 
1 )® = +1 


b) If 7 is odd (that is 7 = 1): 

1. Stays the same as in the previous case. 

2 . Gets an opposite sign, therefore, will cancel out with the a) case 2 in the sum. 


3. Stays the same as in the previous case. 


4. Stays the same as in the previous case. 
Therefore the result is: 


(Ge) = ±^ 2 " 


N—m—1 

E 

^= 0 . 4 .. 



N — m — 1\ 

, /3 + 2 J_ 


om 


( 122 ) 


It is time to fix a sign: if m = 2 mod 8 , the global sign is positive from the Table 6 , and Re[{l + i)^ i] = 2 2 
Therefore, an overall sign is positive. If m = 6 mod 8 , overall sign in negative and Re[{l + 1 )^“™“^] = —2 2 . 

Therefore, an overall sign in positive. □ 


Lemma 8 . A N-qubit (N = Sk A) four-uniform HG state even after tracing out one party violates the classical 
bound by an amount that grows exponentially with number of qubits. Moreover, the violation only decreases with the 
constant factor. 


Proof. Denote M = N — 1. Then 

M 

(Sjv)q = ^ 


m—2,4 

1 

4\/2 


/ \ \M —771+3 


2V2 _ 




1 + 


1 sM / 1 sM 


V2J 


V 2 J 


771 even 

1 nM+3 

vyfJ 


-( 


/ 1 xM+3 


( 123 ) 
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Checking the ratio of the quantum and classical values, we have that 


/ sN-l / sN-l 

{Bn-i)q n^oo + ^ + _ 1.20711^ 

{Bn-i)c 2 ^ ~ ^ 2^2 + 2 ' 


(124) 


While the same value for the = 8fc + 4 qubit four-uniform complete HG state is ^7t°°fa 1 . 20 'ni ^ Therefore, 

after tracing out a single qubit, the local realism violation decreases with the small constant factor. □ 


It is important to note that the similar violation is maintained after tracing out more than one qubit. For example, 
numerical evidence confirms for N = 12, that if one takes a Bell inequality with the odd number of X measurements, 
instead of the even ones as we have chosen in the proof, an exponential violation is maintained after tracing out 2 
qubits. But even if 5 qubit is traced out, the state is stilled entangled and this can be verified using the separability 
inequality [26]. Exact violation is given in the following table: 


#k 

Quantum Value 

Classical Bound 

Separability Bound 

Ri Ratio 

0 

153.141 

64 


2.39283 

1 

89.7188 

32 


2.78125 

2 

37.1563 

32 


1.16113 


3 

15.4219 

- 

72 

10.9049 

4 

6.375 

- 

^J2 

4.50781 

5 

2.70313 

- 

^J2 

1.9114 


Table 7. Violation of Bell (for odd m) and Separability inequalities in iV = 12 qubit 4-uniform HG state. Here k is 
the number of traced out qubits. Red line represents that when k = 1, or equivalently one qubit is traced out, the 
violation of Bell inequalities increases. This is caused by decrease in the classical bound [3]. 


APPENDIX G: SEPARABILITY INEQUALITY VIOLATIONS FOR FULLY-CONNECTED 
THREE-UNIFORM HYPERGRAPH STATES AFTER LOOSING ONE QUBIT 


Lemma 9. The following statements holds three-uniform complete hypergraph states: 
(i) For V = -I- 5, V = -I- 6, or N = ik -\-l: 


(G7) = {X ...XZ ...Zl) = 


m N—m—1 


- 1 ^-1 1 if {m — 1 ) = 0 mod 4, 

-F I I if {m — 1) = 2 mod 4. 

2L~2~ J 


(a) For N = 8k 1, N = 8k 2, or N = 8k -\-3: 


(G7) = {X ...XZ ...Zl) = 


m N—m—l 


-F I I if (m — 1) = 0 mod 4, 




if {m — 1 ) = 2 mod 4. 


(Hi) For N = Ik: 


(Gr) = = 0 

m N—m—1 


(125) 


(126) 


(127) 


Proof. As a starting point, we derive the remaining gates on the right hand side of the equations (125), the same 
derivation turns out to be working for Eq. (126) and (127). This approach is analogous to the Appendix B, but now, 
m is odd. 
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1 . #®® 
2 . #@ 



is odd => C'@@ remains. 


_ (m-l)(m-2) . 
— 2 


even, if {m — 1) = 0 mod 4 
odd, if (m — 1) = 2 mod 4 


C® cancels. 
=> C@ remains. 


3- #{} 

4. # A A 

5. ^ @ A 

6. #A 


_ m(m —l)(m —2) 


2-3 


IS 


even, if (m — 1) = 0 mod 4 
odd, if (m — 1) = 2 mod 4 


gives a positive sign, 
gives a negative sign. 


(™) is odd => Caa remains. 

is even C^a cancels. 


(”) = 


m(m —1) 


2-3 


IS 


even, if (m — 1) = 0 mod 4 => Aa cancels, 

odd, if (m — 1) = 2 mod 4 => Aa remains. 


Table 8 . Counting phase gates for a three-uniform HG when m {m is odd) systems are measured in X direction. 


Consider two cases: 

1. If (m — 1) = 0 mod 4: 


(Gr) = ±{+\^^X_^^2CZ_^l A@@AaaAa«|+)®" 

m N-m-1 V®,A,<C> 

= ±^rr A@®AaaAaoAa ■ 


(128) 


Here ® again refers to X operator, A to Z and <C> to 1 and is denoted by 7 . The strategy is similar to the previous 
case: count the number of -fl’s and -I’s on the diagonal. Their difference divided by 2^, gives the trace. 

We use (—1)'* to define the sign of the diagonal element and s = ( 2 ) + ( 2 ) + ( 1 ) (T) + (i)- ^ even, the value 

on the diagonal is 4-1 and —1, otherwise. We consider all possible values of a, /3 and 7 : 

a) If 7 is even (that is 7 = 0 ): 

Considering the terms from even /3: 

1. /? is even & a is even if a = 0 mod 4 and /3 = 0 mod 4 => (—1)® = 4-1 

if a = 0 mod 4 and (3 = 2 mod 4 => (—1)® = —1 

if a = 2 mod 4 and (3 = 0 mod 4 => (—1)® = —1 

if a = 2 mod 4 and (3 = 2 mod 4 => (—1)* = 4-1 

2. (3 is even & a is odd if (a — 1) = 0 mod 4 and (3 = 0 mod 4 (—1)® = 4-1 

if (a — 1) = 0 mod 4 and (3 = 2 mod 4 (—1)* = —1 

if (a — 1) = 2 mod 4 and (3 = 0 mod 4 (—1)® = —1 

if (a — 1) = 2 mod 4 and (3 = 2 mod 4 (—1)® = 4-1 


Considering odd (3: 
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3. /3 is odd & a is even if a = 0 mod 4 and (/? — 1) = 0 mod 4 => (—1)® = — 1 

if a = 0 mod 4 and (/3 — 1) = 2 mod 4 (—1)® = +1 

if a = 2 mod 4 and (/? — 1) = 0 mod 4 (—1)® = +1 

if a = 2 mod 4 and (/3 — 1) = 2 mod 4 (—1)® = —1 

4. /3 is odd & a is odd if (a — 1) = 0 mod 4 and (/3 — 1) = 0 mod 4 => (—1)® = —1 

if (a — 1) = 0 mod 4 and (/3 — 1) = 2 mod 4 => (—1)® = +1 

if (a — 1) = 2 mod 4 and (/3 — 1) = 0 mod 4 => (—1)® = +1 

if (a — 1) = 2 mod 4 and (/3 — 1) = 2 mod 4 => (—1)® = —1 

b) If 7 is odd (that is 7 = 1): 

Considering the terms from even /3: 

1. and 2. Nothing changes in comparison to a) 1. and 2. 

3. and 4. These two terms have opposite sign from a) 3. and 4. Therefore, in the sum they cancel (7 is always 1.) 
Therefore, 


(G7) = 


N—m — 1 

E 


2^-1 [ V /3 

/3=0,4.. ^ ^ 


N — m — f N — m — \ 


E 

<y—0,4:.. 


+ 


m 


m 


Q. \ J \Q;-t-2/ \Q-t-3 


m 


(129) 


2N- 


^i?e[(l + . (i?e[(l + ir\ + /m[(l + *)-]) . 


Now we can consider each cases separately, for this we use the lookup table 0: 

[i) If TV = 8 A: + 5 and if (m — 1) = 4 mod 8 , Re[{\ + = —2 2 and Re[{l + 1 )"*] + /m[(l + 1 )"*] = 2~^' 

Therefore, 


(Gt) = - 


N-m-2 m + 1 

2 — 2 — . 2— 
2^ 


■© 


(130) 


And, if (m — 1) = 0 mod 8 , Re[{l + i)^ = -\-2 2 and i?e[(l + + /m[(l + i)'^] = —2 2 . Therefore, 


(Gt) = - 


_ N-m-2 _ m + 1 

2 — 2 — . 2^- 


2 N -1 




(131) 


Exactly the same is true for N = 8k + 7. But for iV = 8 fc + 6 and if (m — 1) =4 mod 8 , i?e[(l + i)'^ ™ i] = —2 2 

and i?e[(l + i)"*] + /to[( 1 + i)™] = 2"^^. Therefore, 


( G7 ) = - 


2 ^^. 2 ^ 1 ,^j 


2^-1 




(132) 


And, if (m — 1) = 0 mod 8, i?e[(l + i)^ ™ ^] = 2 ^ ™ ^ and i?e[(l + i)™] + /m[(l + i)™] = —2 "G". Therefore, 


( G7 ) = - 


jv-m-l m + 1 I W-1 I 

2 2 -2 2 /1\L-2-J 


2iv-i 


())' 


(133) 


(iz) If iV = 8 fc + l and if (m—1) = 4 mod 8 , i?e[(l+*)^ ™ ^] =—2 2 and i?e[(l + z)’"]+/m[(l + i)'"] = —2 ^ . 

Therefore, 


( G7 ) 


2^-1 



N-1 

2 


(134) 
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And, if (m — 1) = 0 mod 8, i?e[(l + i) 


N—m—1 


iV “ TTX “ i. 7TL i 

] = 2 2 and Re[{l + i)™] + /m[(l + i)™] = 2^“. Therefore, 


(G7) = 




2N-1 


©' 


(135) 


It it analogous for other two cases as well. 

[iii) For N = 4fc, i?e[(l + = 0. Therefore, {G^) = 0. 


2. If (m — 1) = 2 mod 4: 


(Gr) = C'@@C@CaaCa«CaC^|+)®'^ 

m N-m-1 V®,A,^ 


= —Tr 
2 ^ 


(^00 C(g Ca A Ca<o> C'<o> 


(136) 


Here (S again refers to X operator, A to Z and <0 to 1 and is denoted by 7 . The strategy is similar to the previous 

case: count the number of +l’s and -I’s on the diagonal and their difference divided by 2 ^, gives the trace. 

We use (—1)'* to define the sign of the diagonal element and s = ( 2 ) + ( 1 ) + ( 2 ) + ( 1 ) (T) + (i)- ^ even, the 

value on the diagonal is +1 and —1, otherwise. We consider all possible values of a, /3 and 7 : 

a) If 7 is even (that is 7 = 0): 

Considering the terms from even /3: 

1. /? is even & a is even if a = 0 mod 4 and ,5 = 0 mod 4 => (—1)® = +1 

if a = 0 mod 4 and /3 = 2 mod 4 => (—1)® = —1 

if a = 2 mod 4 and /3 = 0 mod 4 (—1)® = —1 

if a = 2 mod 4 and /3 = 2 mod 4 (—1)® = +1 

2. /3 is even & a is odd if (a — 1) = 0 mod 4 and /3 = 0 mod 4 (—1)® = —1 

if (a — 1) = 0 mod 4 and /3 = 2 mod 4 (—1)® = +1 

if (a — 1) = 2 mod 4 and /3 = 0 mod 4 => (—1)® = +1 

if (a — 1) = 2 mod 4 and (3 = 2 mod 4 => (—1)® = —1 

Considering odd (3: 

3. (3 is odd & a is even if a = 0 mod 4 and (/3 — 1) = 0 mod 4 => (—1)® = +1 

if a = 0 mod 4 and (/3 — 1) = 2 mod 4 => (—1)® = —1 
if a = 2 mod 4 and (/3 — 1) = 0 mod 4 (—1)® = —1 

if a = 2 mod 4 and (/? — 1) = 2 mod 4 ^ (—1)® = +1 

4. (3 is odd & a is odd if (a — 1) = 0 mod 4 and (/? — 1) = 0 mod 4 => (—1)® = —1 

if (a — 1) = 0 mod 4 and {(3 — 1) =2 mod 4 => (—1)® = +1 
if (a — 1) = 2 mod 4 and (/3 — 1) = 0 mod 4 (—1)® = +1 

if (a — 1) = 2 mod 4 and (/? — 1) = 2 mod 4 (—1)® = —1 

b) If 7 is odd (that is 7 = 1):: 

1. and 2. These two terms have opposite sign from a) 1. and 2. Therefore, in the sum they cancel (7 is always 1.) 
3. and 4. Nothing changes in comparison to a) 3. and 4. 
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Therefore, 


(G7) = 


2N-i 


N—m — 1 

E 

/3=0,4.. 


N — m—1 
13 + 1 


N — m—1 
0 + 3 


E 

-I *- (x—0,A.. 


m 

a + 1 


m 

a + 2 


m 

G! H- 3 


(137) 


= ■ (i?e[(l + zr] - /m[(l + *)-]). 

Now we can consider each cases separately, for this we use the lookup table 0: 


(z) If = 8fc + 5 and if (m— 1 ) = 4 mod 8,/m[(l + 7)'^ i] = 2 2 and i?e[(l +7)"*] —/to[(1 + z)™] = —2 ^ . 

And the overall sign from the Table 8 is negative. Therefore, 


(G7) = 


N-m-2 

2—2— 


2N-1 



(138) 


And, if (m — 1) = 0 mod 8, /m[(l + i)^ ^] = —2^ ^ ^ and Re[{l + z)"*] — /m[(l + z)™] = 2”^^^. Overall sign in 

negative. Therefore, 


(G7) 


- 2 " 


2ZV-1 



(139) 


The same holds for TV = 8fc + 6 and N = 8A: + 7. Besides, (zz) differs with the sign flip and is trivial to check.lt is also 
trivial to prove {Hi), as when N = 4fc, Jm[(l + =0. □ 

Lemma 10. A N-qubit (except when N = Ak) three-uniform HG state violates the separability inequality exponentially 
after tracing out a single qubit. 

Proof. We consider only one case, TV = 8fc + 5, as others are analogous. Here M = N — 1: 


m odd ^ ^ 

Separability bound is [26] and it does not depend on the number of qubits. □ 

It is important to note that the similar violation is maintained after tracing out more than one qubit. Numerical 
evidence for TV = 11 is presented below. 


#k 

Quantum Value 

Separability Bound 

Ratio 

0 

511.5 

V2 

361.69 

1 

16 

V2 

11.31 

2 

8 

V2 

5.66 

3 

4 

V2 

2.83 

4 

2 

V2 

1.414 


Table 9. Violation of the separability inequalities [26] in TV = 11-qubit 3-uniform HG states. Here k is the number 
of traced out qubits. When fc = 0, the Mermin inequality is violated as expected. 


[1] P. Hayden, D.W. Leung, and A. Winter, Comm. Math. 
Phys. 265, 95 (2006); D. Gross, S. Flammia, and J. Eis- 
ert, Phys. Rev. Lett. 102, 190501 (2009); M. J. Bremner, 
C. Mora, and A. Winter, Phys. Rev. Lett. 102, 190502 
(2009). 

[2] M. Hein, W. Diir, J. Eisert, R. Raussendorf, M. Van den 


Nest, and H.-J. Briegel, Entanglement in Graph States 
and its Applications, in Quantum Computers, Algorithms 
and Chaos, edited by G. Gasati, D.L. Shepelyansky, P. 
Zoller, and G. Benenti (lOS Press, Amsterdam, 2006), 
quant-ph/0602096. 

[3] N. D. Mermin, Phys. Rev. Lett. 65, 1838 (1990). 

[4] O. Giihne, G. Toth, P. Hyllus, and H. J. Briegel, Phys. 
Rev. Lett. 95, 120405 (2005); G. Toth, O. Giihne, H. J. 
Briegel, Phys. Rev. A 73, 022303 (2006). 

























31 


[5] O. Giihne and G. Toth, Phys. Rep. 474, 1 (2009). 

[6] C. Kmszynska and B. Kraus, Phys. Rev. A 79, 052304 
(2009). 

[7] R. Qu, J. Wang, Z. Li, and Y. Bao, Phys. Rev. A 87, 
022311 (2013). 

[8] M. Rossi, M. Huber, D. Brufi, and G. Macchiavello, New 
J. Phys. 15, 113022 (2013). 

[9] O. Giihne, M. Guquet, F. E. S. Steinhoff, T. Moroder, M. 
Rossi, D. Brufi, B. Kraus, and C. Macchiavello, J. Phys. 
A: Math. Theor. 47, 335303 (2014). 

[101 X.-Y. Chen and L. Wang, J. Phys. A: Math. Theor. 47, 
415304 (2014). 

[11] D. W. Lyons, D. J. Upchurch, S. N. Walck, and C. D. 
Yetter, J. Phys. A: Math. Theor. 48, 095301 (2015). 

[12] M. Rossi, D. Brufi, and C. Macchiavello, Phys. Scr. 
T160, 014036 (2014). 

[13] H. Buhrman, R. Cleve, J. Watrous, and R. de Wolf, Phys. 
Rev. Lett. 87, 167902 (2001), C.E. Mora, H.J. Briegel, 
and B. Kraus, Int. J. Quantum Inform. 5, 729 (2007). 

[14] A. B. Grilo, 1. Kerenidis, and J. Sikora, arXiv: 1410.2882. 

[15] B. Yoshida, arXiv: 1508.03468. 

[16] J. Miller and A Miyake, arXiv: 1508.02695. 

[17] L. Hardy, Phys. Rev. Lett. 71, 1665 (1993). 

[18] G. Svetlichny, Phys. Rev. D 35, 3066 (1987). 

[19] S. Popescu and D. Rohrlich, Phys. Lett. A 166, 293 
(1992). 


[20] R. F. Werner and M. M. Wolf, Phys. Rev. A 64, 032112 
(2001); M. Zukowski and C. Brukner, Phys. Rev. Lett. 
88, 210401 (2002). 

[21] N. Brunner, D. Cavalcanti, S. Pironio, V. Scarani, and 
S. Wehner, Rev. Mod. Phys. 86, 419 (2014). 

[22] C. Brukner, M. Zukowski, J.-W. Pan, and A. Zeilinger, 
Phys. Rev. Lett. 92, 127901 (2004). 

[23] S. Abramsky and C. Constantin, Electronic Proc. Theor. 
Comp. Science 171, 10 (2014), arXiv:1412.5213. 

[24] Z. Wang and D. Markham, Phys. Rev. Lett. 108, 210407 

( 2012 ). 

[25] The appendix can be found in the supplementary online 
material. 

[26] S. M. Roy, Phys. Rev. Lett. 94, 010402 (2005). 

[27] V. Giovannetti, S. Lloyd, and L. Maccone, Science 306, 
1330 (2004). 

[28] W.-B. Gao, C.-Y. Lu, X.-C. Yao, P. Xu, O. Giihne, A. 
Goebel, Y.-A. Chen, C.-Z. Peng, Z.-B. Chen, and J.-W. 
Pan, Nature Phys. 6, 331 (2010). 

[29] M. J. Hoban, E. T. Campbell, K. Loukopoulos, and D. 
E. Browne, New J. Phys. 13, 023014 (2011). 

[30] D. Collins, N. Gisin, S. Popescu, D. Roberts, and V. 
Scarani, Phys. Rev. Lett. 88, 170405 (2002). 

[31] M. Koashi ,V. Buzek, and N. Imoto, Phys. Rev. A 62, 
050302(R) (2000) 


